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Abstract 

Data collected with the BABAR detector at the PEP-II asymmetric B Factory at SLAC are used 
to study the lifetimes of the B'^ and B~^ mesons. The data sample consists of 7.4 fb~^ collected 
near the T(45) resonance. B^ and B~^ mesons are fully reconstructed in several exclusive hadronic 
decay modes to charm and charmonium final states. The B lifetimes are determined from the 
flight length difference between the two B mesons which are pair-produced in the T(4S) decay. 
The preliminary measurements of the lifetimes are 

r^o = 1.506 ± 0.052 (stat) ± 0.029 (syst) ps, 
tb+ = 1.602 ± 0.049 (stat) ± 0.035 (syst) ps 

and of their ratio is 

Tb+/tbo = 1.065 ± 0.044 (stat) ± 0.021 (syst). 
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1 Introduction 



Precise measurements of the charged and neutral B meson hfetimes are needed for tests of theo- 
retical models of heavy quark decays. An uncertainty of around 1% will be required to distinguish 
between models |2| predicting different lifetime hierarchies among the B hadrons. The life- 
times are also required for other measurements and they provide a check of our understanding of 
the detector. The current world averages Q for the B lifetimes are r^o = 1.548 it 0.032 ps and 
Tb+ = 1.653 ± 0.028 ps and for their ratio is Tb+/tbo = 1.062 ± 0.029. 

This paper presents preliminary measurements of the B^ and B^ /B^ lifetimes and their ratio 
performed with data collected by the BABAR detector at the PEP-II asymmetric B Factory. At the 
T(4S'), B mesons are produced in B^B^ and B^B^ pairs. In this analysis, one of the B mesons is 
fully reconstructed in a variety of clean hadronic two-body decay modes to charm and charmonium 
final states. An inclusive technique is used to reconstruct the decay vertex of the second B in an 
event, and the lifetimes are determined from the distance between the decay vertices of the two 
B mesons. This novel method, developed for use at an asymmetric B Factory, deals with event 
topologies similar to those addressed in the analyses of CP asymmetries where the time dependence 
is obtained from the distance between the two vertices. The measurement of B lifetimes using this 
technique is therefore a validation of the CP analyses. The systematic errors are quite different 
from previous measurements of the lifetimes and we expect substantial increases in the statistical 
precision in the future. 

The paper is organized as follows. The method used to determine the B meson lifetimes 
is outlined in Section |2|. The data sample and the detector components most relevant for this 
measurement are reviewed in Section ^ The event selection criteria and sample composition are 
described next in Section ^. The determination of the distance between the decay vertices of the two 
B mesons in a selected event is discussed in Section ^ together with the corresponding resolution 
function. The fitting procedure used to extract the B lifetimes from the neutral and charged 
B meson samples is described in Section |6|. Systematic uncertainties are discussed in Section |^. 
The results are given at the end of the paper. 



2 The decay length difference technique 

The T(45) decays exclusively into pairs of charged or neutral B mesons. At the PEP-II asymmetric 
B Factory the T(45) decay products are Lorentz boosted {P^ ^ 0.56)|^ and travel far enough for 
their flight paths (~ 260 /xm) to be comparable to the experimental resolution (~ 150 /im). Since no 
charged stable particles emerge from the T(4S') decay point, the production point of the B mesons 
is unknown and the B lifetimes have to be determined solely from their decay points. The difference 
between the decay lengths of the two B mesons in an event has been used for that purpose in this 
analysis. 

Throughout the paper the signed quantity Az = z^^c — -Zopp is used for the projection on the 
z axis (defined below) of the decay length difference between the fully reconstructed B meson, 
thereafter called -Breo and the opposed B meson, -Bopp- The event topology is sketched in Fig. |l|, 
which is not drawn to scale. In the simplified situation where the B mesons are produced exactly 
at threshold and the T{AS) moves exactly along the z axis, |Az| is distributed exponentially with 
an average of (|Az[) = {I3^)bctb = {PT{AS)/^T(iS))(^^B- Complications arise from several effects. 
The axis of the PEP-II beams is tilted by 20 mr with respect to the z axis, which is defined as 

^The average boost of the center of mass frame is measured on a run-by-run basis using two- prong events M]. 
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parallel to the solenoidal magnetic field, and the energies of the beams fluctuate, giving the T{4S) 
momentum a Gaussian distribution with a standard deviation of 6 MeV/ c. Furthermore, the energy 
release in the T(4S) BB decay makes the B mesons move in the T{AS) rest frame, resulting in 
a non- vanishing opening angle (smaller than 214 mr [^) between the trajectories of the B mesons. 
The latter effect has the highest impact on the Az distribution. In practice, all of these effects are 
small compared to the experimental resolution on Az. 




Figure 1: Event topology. The figure is not drawn to scale (see text). 



3 The BaBAR detector and data set 

The data used in this analysis were collected by the BABML detector at the PEP-II storage ring 
in the period from January to June, 2000. The total integrated luminosity of the data set is 
7.4 fb"^ collected near the T(45) resonance and 0.9 fb~^ collected 40 MeV below the resonance. 
The corresponding number of produced BB pairs is estimated to be 8.4 x 10^. 

The BABAR detector is described elsewhere [Q]. For this measurement, the most important 
subdetectors are the silicon vertex tracker (SVT) and the central drift chamber (DCH), which jointly 
provide charged particle tracking, and the Csl electromagnetic calorimeter (EMC), from which 
photons and vr'^s are reconstructed. The charged particle transverse momentum {pt) resolution is 
approximately described by the formula {5pt/pt)^ = (O.OOISpt)^ + (0.005)^, where is in GeV/c. 
The SVT, with a typical single-hit resolution of 10 //m, provides vertex information in both the 
transverse plane and in z. The precision on charged particle momenta, neutral particle energies and 
spatial coordinates from the tracking and calorimetry leads to resolutions on invariant masses and 
other kinematical quantities which are adequate for a clean separation of the exclusive hadronic 
modes as described in Section ^ Impact parameter resolutions in the transverse and longitudinal 
coordinates are ~ 50 /im at high momentum, and better than 100 //m for px > 0.6GeV/c. The 
achieved precision on exclusive and inclusive vertices are given in Section |5[ 
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Particle identification is needed for complex decay modes and J/il) reconstruction. Leptons 
and hadrons are identified using a combination of measurements from all BABAR components, 
including the dB/dx energy loss using a truncated mean of 40 samples (maximum) in the DCH 
and 5 samples in the SVT. Electrons and photons are identified in the barrel and forward regions 
by the EMC. Muons are identified in the instrumented flux return (IFR). In the central polar 
angle region the Cherenkov ring imaging detector (DIRC) provides a separation better than three 
standard deviations between pions and kaons over the upper part of the momentum range of B 
decay products (above 700MeV/c) where the dB/dx sensitivity vanishes. 

4 Event selection 

B^ and B^ mesons are reconstructed in the following hadronic modesQ; B'^ — > L'(*)^7r+, D^*^^ p'^ , 

D^*'>^at, J/i)K*^ and B+ 1)^*^%+, J/^K+, ^{2S)K+ . Ah final state particles are recon- 
structed. The selection criteria are devised to obtain final samples with approximately 90% purity 
and have been studied using the Monte Carlo simulation. Optimization for background rejection 
was performed using continuum data and sideband regions that do not contain signal events. The 
details of the selection of these events are described below. 

4.1 Selection criteria 

The basic objects used in the selection are charged tracks and neutral electromagnetic clusters. 
Charged tracks are built out of hits in the SVT and the DCH by a Kalman filter-based tracking 
algorithm Q. To be retained, a track must not miss the average beamspot position by more 
than 1.5 cm in xy and 10 cm in z. Most tracks must penetrate halfway radially into the drift 
chamber (at least 20 hits are required) and have a momentum between 0.1 and lOGeV/c. The 
latter restrictions do not apply to soft pions or low momentum tracks from a high multiplicity 
D decay. Electromagnetic clusters in the EMC are photon candidates when they carry more than 
30 MeV, have a narrow transverse shower profile (the lateral moment defined in Ref. @) less than 0.8 
and are not matched to a charged track. 

The selection procedure begins by enforcing standard criteria 0] on the charged multiplicity, 
total energy, global event vertex position and event shape which enrich the sample with multi- 
hadronic events from T(4S') decays. A tight requirement on the normalized second Fox- Wolfram 
moment i?2(= H2/Hq) < 0.5 is applied in all analyses pertaining to this paper to reject jet-like 
events from the continuum. 

The selection then proceeds by building a tree of composite objects (fast decaying particles) 
from tracks and photon candidates, culminating in a. B candidate if one exists in the event. Geo- 
metrical fits are used to reconstruct decay points and declared successful when the probability 
is above 0.1%. No goodness of fit requirement is enforced for the kinematic fits used to impose 
mass constraints. 

Neutral pion candidates are formed from pairs of photon candidates with an invariant mass 
within 20 MeV/c^ of the nominal tt^ mass and with a minimum energy of 200 MeV, except for those 
used in attempts to reconstruct a D*^. A kinematic fit is performed on the selected candidates 
imposing the vr*^ mass constraint. 

Kg candidates are reconstructed in the two charged pion decay mode only. Oppositely charged 
tracks are vertexed using a geometrical fit. The tt+tt^ invariant mass, computed at the vertex of the 

^Throughout this paper, when a mode is quoted, its charge conjugate is imphed. 
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two tracks, has to lie between 462 and 534MeV/c^. The opening angle between the flight direction 
and the candidate momentum vector must be smaller than 200 mr. Finally, the transverse 
flight distance from the event primary vertex is required to be greater than 2 mm. 

K and vr meson candidates are combined to search for the following D meson decay modes; 
D° K-7r+, /C-vr+vrO, K-tt+tt+tt-, K^tt+tt- and D+ K-tt+tt+, K^^tt+. All kaons and the 
pions from two body modes are required to have a momentum above 200MeV/c. This criterion is 
relaxed to 150MeV/c for pions from three or four body final states. Candidates with an invariant 
mass within three standard deviations of the known mass are retained. The standard deviation 
computed event by event by the tracking algorithm was used after it was verified that the normalized 
error distribution (pull) was correctly estimated at the 10-20% level. The reconstruction of — >■ 
K^TT^TT^ is restricted to the dominant resonant mode K~p~^, p+ tt+tt^. The n'^ir^ mass 

is required to lie within ibl50MeV/c^ of the known p mass and the p decay helicity angle in its rest 
frame, 0*, is required to satisfy |cos0*| > 0.4. All D candidates must have a momentum greater 
than 1.3 GeV/c in the T(45') frame and have a successful vertex fit. The accepted D candidates are 
mass constrained before they are used in subsequent intermediate state searches. Charged (neutral) 
D* mesons are built by combining a with a soft charged (neutral) pion with momentum above 
70 (100) MeV/c. Constraining the vertex of the D*~^ to the beamspot vertical position within 40 pm 
improves the resolution on the soft pion direction and hence on Am, the difference between the D* 
and the D candidate invariant masses. This procedure has been verified to be bias-free in Monte 
Carlo events. An acceptable range of ±3 standard deviations around the nominal Am completes 
the D* selection. The range limits (zbl.lMeV/c^ for K' p+ , zbO.SMeV/c^ for ah other D*+ 

modes and ±1.4MeV/c2 for the D*^) are weighted averages of the standard deviations obtained 
from two-Gaussian fits to the experimental Am distributions. 

The charmonium meson selection is the same as that described in Ref. Q. Leptonic decays of 
charmonium mesons are reconstructed from pairs of oppositely charged tracks considered either as 
electron or muon pairs. At least one decay product must be positively identified as a lepton. For 
an electron inside the EMC acceptance, this means a positive electromagnetic shower signature, 
otherwise that the measured dE'/dx be consistent with the electron hypothesis. For muon pairs, 
the second track has to be minimum ionizing in the EMC. In the electron case, when both tracks 
are inside the EMC acceptance and are such that E/p > 0.5, where E is the energy measured in 
the EMC and p is the momentum obtained from the tracking, bremsstrahlung photons found in the 
EMC are combined with the tracks to reconstruct the electrons prior to radiation. J/tp e~^e~ 
{p'^p~) candidates must have an invariant mass between 2.95 (3.06) and 3.14MeV/c2. The invariant 
mass of the '4^{2S) — > p'^p~ candidates is required to be within 50MeV/c^ of the known 1^(23) 
mass, but a less stringent lower limit of 200MeV/c^ is imposed for 'ip{2S) — > e'^e~ candidates. 
ip{2S) — > J/'ip-7T~^TT~ decays are also reconstructed by combining oppositely charged pions with mass- 
constrained J/^p candidates and selecting the ip{2S) invariant mass between 3.1 and 4.0GeV/c^ and 
the tp{2S) — J/iIj mass difference in a window of ±15MeV/c^ around the nominal value. 

B candidates are formed by combining a D* , D, J/ip or ip{2S) candidate with a vr"*", /j"*", 
af, K*^ or candidate having a minimum momentum of 500MeV/c in the T{4S) frame. For 
B'^ D^*^^ p'^ , the TT^ meson from the p"*" decay is required to have an energy higher than 300 MeV. 
For B^ — D^*^~af, the af meson, selected by combining three charged pions, must have an 
invariant mass between 1.0 and 1.6 MeV/ and a successful vertex fit. 

Kaon identification is used to reject background at the price of some loss in reconstruction 
efficiency. For most of the B^ modes, a loose kaon selection, or no selection, is enough to reduce 
the background level to the required level. Tighter kaon identification is mandatory only for modes 
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with high combinatorial background, such as i? D^a^ . 

Finally, to achieve the required high signal to background ratio, stricter requirements on event 
shape variables are applied. A selection is performed on the angle, 9t, between two thrust axes, 
namely the thrust axis computed from the tracks of the B meson candidate and the thrust axis of 
all the other charged tracks and neutral clusters in the event. A flat distribution is expected for 
B meson pairs while continuum events tend to peak at cosOt = ±1. For the B — > D*p modes, 
I cos6t\ is required to be less than 0.9 for the K^tt^ and K~tt^tt^ channels and less than 0.8 

for K~TT~^TT~^TT~ and K'^tt'^tt^ . The requirements are 0.8 and 0.7 for the B^ — > D*~af modes 

and the same two couples of D modes. There is no restriction based on thrust for the B — > D*~n~^ 
decays. For the B^ decaying to charged D mesons, the requirements are 0.9, 0.8 and 0.7, if the D~ 
is accompanied by a 7r+, a p+ or an af. For the B^ D^*^^7r^ decays, | cos 6*^1 is required to be 
less than 0.9 in the case of the K^tt^ decay and to be less than 0.8 in all other cases. The 

requirement is 0.9 for the B decays to charmonium. 

4.2 Sample composition 

Signals for each decay mode in the selected sample are isolated by analyzing the two-dimensional 
distribution of the kinematical variables, AE and mssi which have an intuitive meaning in the 
T(45) rest frame. AE = E*^^ — E^ is the difference of the B candidate energy and the beam energy. 

ruES = \J ~ Prec the mass of a particle with a reconstructed momentum p*^^ = J^i Pi assumed 
to have the beam energy, as is the case for a true B meson. Candidates with 5.2 < m-Es < 5.3 MeV/ 
and \ AE\ < 140 MeV are considered in further analyses. In case an event has several B candidates, 
only the one with the smallest AE' is kept. The AE' and m-Es variables are essentially uncorrelated. 
The resolution on m-Es is about 3MeV/c^ and is dominated by the energy spread of the beams, 
inherent to the collider. The resolution on AE varies from mode to mode between 12 and 40 MeV. 
For each mode a rectangular signal region is defined by the three standard deviation bands in rriES 
(5.27 < rriEs < 5.29GeV/c^) and AE (an interval which depends on the mode). 

For each mode, the sample composition was determined by fitting the tit-es distribution, for 
candidates within the signal region in AE, to the sum of a single Gaussian representing the signal 
and a background function used by the ARGUS collaboration |^. The purity of each subsample 
is computed as the ratio of the area of the Gaussian in the ±3a range over the total area in this 
range. The results are summarized in Table |l[ Figure Q shows the niEs distributions for the 
summed hadronic B^ and B~^ modes with the fits superimposed. 

To extract separately the B^ and B~^ lifetimes, a thorough understanding of the background 
shape and composition is necessary. In the signal region, the background receives contributions 
from continuum, B^B^ and B^B^ pairs in various proportions for the different modes. A consistent 
background shape is found from studies of several samples; namely, on-resonance data events outside 
the signal region, off-resonance data, and generic BB and continuum Monte Carlo samples. The 
contamination of the B^B~ subsamples by B^B^ events, and vice versa, is found to be small. 

5 Vertex reconstruction and resolution function 

To extract the B meson lifetimes from the sample, the separate decay vertices have to be identified 
and the distance between them along the z axis, Az, evaluated. The resolution function for Az 
is the most critical ingredient of this analysis as the width of the Az distribution is due to the 
combined effect of the B lifetime and the detector resolution. 
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Table 1: Two-body hadronic and B^ decay candidate yields and signal purities from the fit to 
the mES distribution. Signal purities are estimated for mss > 5.27GeV/c^. 



Decay 


mode 


Number of B candidates 


Purity (%) 


B^^ 
B^ ^ 
5° ^ 
B^ ^ 
i?o ^ 
5° ^ 
5° ^ 
5+ ^ 


D*-TT+ 

D*-p+ 

D*-af 

D-t:+ 

D-p+ 

D-af 

DOtt+ 

J/i{j,tp{2S)K+ 


552±26 
374±23 
202±18 
537±25 
279±20 
194±18 
167±15 
1528±43 
446±25 
294±17 


90 
84 
79 
90 
84 
73 
90 
88 
89 
99 


Total 
Total B+ 


2210±58 
2261±53 


86 
89 



A geometric and kinematic fit of the fully reconstructed Srec is performed and is required to 
converge. The masses of D and J/tl' resonances are fixed at the nominal values Furthermore, the 
vertices of short-living resonances are constrained to be identical to the B vertex. Two-Gaussian 
fits adequately describe the error distributions of the full Monte Carlo simulation used to model 
the signal events, with standard deviations ranging from 45 to 65 pm and normalized fiuctuations 
(called pulls in the rest of the paper) around 1.1. 

The vertex of the other B is determined via an inclusive and iterative procedure applied to all 
the tracks that are not associated with the -Brec- To avoid biases from long-lived particles such 
as neutral strange particles decaying to two charged prongs (a topology known as a V" decay), 
as many V'' decays as possible are reconstructed and subsequently treated as single tracks. A 
common vertex is then fitted using all tracks including the -Bopp "pseudotrack" which is obtained 
from the T(4S') decay kinematics using the Srec and the beamspot. The procedure is repeated, after 
removing tracks and V*^ decays which result in poor fits, until stability is achieved. The algorithm is 
highly efficient and successful in minimizing biases due to the secondary charm decays of the -Bopp- 
From the weighed means of two-Gaussian fits to Monte Carlo distributions of the residuals, one 
standard deviation errors of the order of 115 pm and biases around 25 pm are obtained. Although 
still present, the bias due to the charm particle lifetime is kept small. 

With both vertices reconstructed, the decay length difference Az is obtained. The Monte Carlo 
simulation was used to study the Az resolution function 7^ as a function of two variables, the residual 

generated and the pull 5( Az)/(j( Az) . Figure y presents the distributions for 
the decay B~^ J/ipK'^. Two-Gaussian fits successfully describe the Monte Carlo residual and 
pull distributions. They give a one standard deviation width of 130 ± 1 pm and a bias of 24.5 =b 1 ^m 
for the residual. The corresponding values for the pull are 1.21 it 0.01 and 0.29 ± 0.01 respectively. 
As expected their widths are dominated by the i?opp vertex errors. Tails remain due to charm 
decays of the -Bopp- Because it is governed by the -Bopp vertex, the resolution function shape is the 
same for all modes and also for neutral and charged B mesons. These comparisons are summarized 
in Table ^. 
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Figure 2: mES distribution for all the hadronic modes for (a) and (b) B^. The signal is fitted 
to a Gaussian and for background we use the ARGUS parameterization 0]. The total numbers of 
signal events in all B^ and B^ modes are 2210 it 58 and 2261 it 53 respectively. 

A common TZ function can be used for all modes. The shape of the resolution function can be fit 
with several parametrizations. Results for the pull representation (which is used in the likelihood 
function of the lifetime fit) are summarized in Table |^. It can be seen that a single Gaussian is too 
simple whereas a two-Gaussian line shape G + G fits adequately. A fair agreement is also obtained 
with a function of the form G + G E, the sum of an unbiased Gaussian G and of the convolution 
of the same Gaussian with a decaying exponential E, as can be seen in Fig. ^. 

Table 2: Parameters of the G + G ® E resolution function TZ in the pull representation for a few 
B^ and B^ decay modes. The parameter g is the fraction of events in the pure Gaussian, s is the 
standard deviation of the Gaussian function and is the decay constant of the exponential. A 
single TZ function for all modes is used in the lifetime fit. 





g S Tr 


5° ^ D*-t:+ 
B+ J/'4)K+ 
B+ ^ ^Ovr+ 


0.637 ±0.031 1.019 ±0.017 0.914 ± 0.062 
0.666 ±0.017 1.047 ±0.010 1.024 ±0.041 
0.678 ±0.015 1.033 ±0.008 0.956 ± 0.035 
0.644 ±0.020 1.008 ±0.010 0.833 ± 0.037 



In a higher statistics Monte Carlo study, tails become apparent which cannot be reproduced by 
any of the parametrizations considered above. These "outlier" events in the tail are bound to bias 
the lifetime measurement if they are not accounted for in the fit. Some are due to long lived strange 
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Table 3: Values of the parameters for the different functional forms used to fit the Az resolution 
function. 



Parametrization 


Parameters 


X^/ndof 


Single Gaussian (G) 


mean = 0.273 ± 0.017 
width = 1.148 ± 0.015 


135.3/51 


Gaussian centered at zero plus 

the same Gaussian convoluted with an 

exponential {G + G ^ E) 


g = 0.643 ± 0.029 
s = 1.012 ± 0.0164 
Tr = 0.936 ± 0.062 


62.8/50 


Two Gaussians [G + G) 


/ = 0.808 ± 0.047 
si = 0.998 ± 0.032 
hi = 0.204 ± 0.025 
S2 = 1.931 ± 0.143 
62 = 0.819 ± 0.139 


40.0/48 




Figure 3: Distribution of (a) residual and (b) pull of Az from a sample of J/ipK'^ Monte 

Carlo simulated events. 
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Figure 4: Fit of the Az pull distribution on a sample of 5000 — > D^n^, — > A'+vr^vr" 
Monte Carlo simulated events, with a function consisting of a Gaussian centered at zero and the 
same Gaussian convoluted with an exponential {G + G ® E parametrization) . 

particles forming a vertex with, e.g., the -Bopp pseudotrack which is robust enough to survive the 
iterative vertex finding algorithm. To eliminate outliers, a requirement of a minimum of two real 
tracks to form the -Bopp vertex, |Az| < 3000 /um and it(Az) < 400 ^m is imposed. There remain 
1752 lb 46 (1879 ± 45) signal events in the neutral (charged) B lifetime samples after the selection. 

The distribution of a{/S.z) can be fitted to a Crystal Ball line shape function The fitted 
function is denoted p{(j). Use of that spectrum with the pull representation of the resolution 
function provides an elegant way of implementing the lifetime fits as explained in Section |6|. 

Figure ^ shows the Az distributions in the data for candidates satisfying tties > 5.27GeV/c^. 
Superimposed are the Monte Carlo simulation predictions using the central values of the lifetimes 
of the B mesons from Ref. Ol. 



6 Lifetime fits 

The B meson lifetimes are extracted from the Az distributions of the selected events with an 
unbinned maximum likelihood fit. The measurements performed on each event i are represented by 
three input numbers: (Az)^, its error, cjj, and Psig,i- All events which satisfy 5.2 < tties < 5.3 GeV/ 
and have IS.E in the signal region are input to the fit. Psig,i is the probability for an event i to come 
from signal. This probability is obtained from the independent fit of the rriEs spectrum using the 
ARGUS background function and a single Gaussian, as described previously (see Fig. |2|). The fit 
determines the B lifetime, tb, the proportion of outliers, Poutj and the two sets of parameters ^resoi 
and ^bkg which describe the resolution function and the background shape, respectively. 

The total probability density function (PDF) for an event i, is the sum of three terms 
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Figure 5: Comparison of At distributions in the data (points) and BB Monte Carlo events (solid 
line histogram) for (a) B^ and (b) B^. The contribution from the other type of -BS-decays in 
each case and from continuum events is also shown (hatched histogram). The simulation has been 
normalized to the data statistics. 

describing the signal, the background and the outliers, 

Ti = p{ai) X [psig,j S{{Az)i,ai;TB,0,csoi) + (1 - Psig,i - Pont) B{{Az)i;9y,i,g) + Pout 0{{Az)i)] . 

The signal PDF, S, is the result of a convolution of the theoretical Az distribution (two expo- 
nential wings) and the pull representation of the Az resolution function. Three or five resolution 
parameters, ^resoh (depending on which parameterization is chosen for TZ, G + GotG + G^E) are 
used to model the width and non-Gaussian tails. The parameters are free in the lifetime fit, i.e. the 
values of the parameters of the resolution function are extracted from the data themselves. Reduc- 
ing the number of parameters modeling TZ, which are unavoidably correlated to tb, has the effect 
of improving the statistical precision on the lifetime. The parameterization needs to be flexible 
enough to reproduce all features of the resolution function, however. The optimal choice depends 
on the statistics available. 

The events in the substituted mass sideband carry most information on the background Az dis- 
tribution. A satisfactory description of this distribution is given by a function adding a single 
Gaussian and two independent exponential tails for the negative and positive Az ranges. This pa- 
rameterization uses six free parameters ^bkg^ the width and mean of the Gaussian, the "lifetimes" of 
the two exponentials and the fraction of events in each of the exponentials. Such a shape is flexible 
enough to account for the combinatorial background with continuum and BB contributions. 

The outlier PDF, O, is a single wide Gaussian with a fixed width of 2500 //m and a fixed mean 
of zero. 

A full Monte Carlo sample with 5000 events in the channel B^ D~it~^ , D~ K~^it~-k~ was 
used to compare different parameterizations of the resolution function (see Table P) which were 





17 



injected into signal-only lifetime fits. The G + G model of the resolution function gives an unbiased 
result but the statistical error on the lifetime is increased by more than 50% compared to the case 
of a fixed TZ. The single Gaussian model is both biased (overestimation of ~ 1.5%) and gives a 30% 
increase of the statistical error. The G + G ® E is intermediate with a 0.3% overestimation of tb 
and the smallest increase of the statistical error (10%). With the present statistics, the G + G ® E 
model of TZ lead to the smallest overall error. With high statistics, the G + G model will become 
the optimum method. 

To determine the B lifetimes, a combined fit to the Az distributions of the and B~^ samples is 
performed. The two Az distributions are not combined but are fitted simultaneously, with different 
sets of parameters (one per charge) to describe the lifetime, the background and the outliers. Since 
the resolution functions for neutral and charged B mesons are compatible, a single set of parameters 
is used to describe the resolution function for both samples. The result of the fit is 

r^o = 1.512 ±0.052 ps, 
tb+ = 1.608 ±0.049 ps, 

where the errors are statistical only. The lifetime fit results are superimposed upon the B^ and 
B^ IS.Z distributions from the B reconstruction signal region in Figs. ^ and ^. The lifetime ratio 
is determined with a similar combined fit, this time to the B^ lifetime and the ratio, which results 
in a value of 

^b+Itbo = 1.065 ±0.044, 

where, again, the error is statistical only. The final results after corrections described in the section 
on systematic errors are applied, are quoted in section ^. 

7 Systematic uncertainties 

The systematic errors associated with event selection, A.z measurement and lifetime fitting are 
discussed in turn. 

It has been verified, using the B^ and B^ Monte Carlo samples, that the event selection does not 
distort the generated Az distribution so that the lifetime measured from these samples is consistent 
with the input value. These checks are done with a limited number of simulated events which results 
in a statistical error on the lifetimes extracted from the fit. This Monte Carlo statistical error of 
1.0% for S° and 0.9% for B+ is taken as a systematic error on the lifetime measurements, labeled 
as "MC statistics" in Table |. 

The Az resolution function 7^ is a key ingredient of the lifetime measurements. The parameters 
describing IZ are results of the fit. Since they are correlated with the lifetime, part of the systematic 
error which comes from their uncertainty is transferred into the statistical error on the lifetime. 
This effect is estimated to be about 2%. For the systematic error associated with the resolution 
function modeling, two contributions are quoted. One is the difference between lifetimes obtained 
with fits using the G + G and G + G ® E parametrizations, a 0.3% effect. The other, at the 
level of 0.3% as well, is an overestimation of the distortion which may result from the simplifying 
assumption that a unique resolution function is suitable for all the modes. 

The sensitivity of the lifetime measurement to quantities which enter into the determination of 
the -Bopp vertex, namely the beam spot position and vertical size, and the -Brec momentum vector 
has been found to be negligible. 
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Figure 6: At distribution for /B^ candidates in the signal region. The result of the lifetime fit 
is superimposed. The background is shown by the hatched distribution. 
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Figure 7: Az distribution for B candidates in the signal region. The result of the lifetime fit is 
superimposed. The background is shown by the hatched distribution. 
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The fit accounts for outlier events via a wide Gaussian function with fixed mean and width. 
Increasing the width from 2500 to 4000 fim gives a change of 1.1% (1.3%) to the {B^) hfetimes. 
A change of 1000 fim in the mean gives an effect smaller by a factor three. 

The effects of the detector geometry and alignment have been studied. These are smaller than 
the effects due to the charm lifetimes, and both parameterizations of the resolution function are 
flexible enough to accommodate small changes in the detector response. The associated systematic 
errors are included in those given for the resolution function model. 

An overestimation of the uncertainty on the length scale measured along z has been obtained 
in considering the mechanical integrity of the SVT. A systematic error of 1% is assumed. Actual 
measurements of the length of the beam pipe with tracks are in progress and should considerably 
reduce this contribution to the error. 

Systematic effects due to the boost enter due to the approximation of using the T(45') boost 
for both B mesons in an event, as well as measurement error on the boost itself. The effect of the 
approximation is estimated to be a 0.4% shift and the measurement error from two-prong events 
is 0.3%. The central values for the lifetimes obtained from the fit are corrected for this shift and a 
conservative error of 0.4% is assigned. Adding the two contributions in quadrature, the systematic 
error due to the boost effects is 0.5%. 

The systematic error due to the signal probability, psig, is obtained by propagating the one 
standard deviation errors on the parameters of the line shape function to the lifetime. A 0.2% 
uncertainty is obtained. 

The uncertainty in modeling the background is estimated by comparing the Monte Carlo dis- 
tributions for the background in the signal region with the mss sideband, and also comparing 
rriES sidebands in data and Monte Carlo simulation. The disagreement is more pronounced for B^ 
mesons (1.1%) than for B^ mesons (0.4%). 

The stability of the results has been checked by splitting the data sample into subsamples of 
different Bj-^c modes and different times of data taking. A summary of all systematic errors as well 
as their sum in quadrature is given in Table ^. 



Table 4: Summary of the systematic errors for the B lifetime measurements. 



systematic effect 


S{tbo) 


5(tb+) 




comment 




(fs) 


(fs) 






MC statistics 


16 


14 


0.014 




Parameterisation of resolution function TZ 


5 


5 


0.001 


some included in the 










statistical error, see text 


One single resolution function 


5 


5 


0.008 




Beamspot, p_b„^ 










Az outliers 


16 


20 


0.005 




Geometry and alignment 








included in TZ 


z scale 


15 


16 




cancels in ratio 


Boost 


6 


6 




cancels in ratio 


Signal probability 


3 


2 


0.005 




Background modeling 


5 


17 


0.011 




Total in quadrature 


29 


35 


0.021 
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8 Summary 



The preliminary results for the B meson lifetimes are 

tbo = 1.506 ± 0.052 (stat) ± 0.029 (syst) ps, 
Tb+ = 1.602 ± 0.049 (stat) ± 0.035 (syst) ps 

and for their ratio is 

tb+/tbo = 1.065 ±0.044 (stat) ±0.021 (syst). 

These results are consistent with previous B lifetime measurements Q and competitive with the 
most precise ones. The statistical errors dominate the uncertainties, while the systematic uncer- 
tainty on the lifetime ratio is smaller than the overall error of the world average for this quantity. 
There are ways to reduce some contributions; however, to reach a 1% total error on the B lifetime 
ratio remains a challenge. 
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